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In theories of spin-dependent radical pair reactions, the time evolution of the radical pair, including 
the effect of the chemical kinetics, is described by a master equation in the Liouville formalism. For 
the description of the chemical kinetics, a number of possible reaction operators have been formu- 
lated in the literature. In this work, we present a framework that allows for a unified description of 
the various proposed mechanisms and the forms of reaction operators for the spin-selective recom- 
bination processes. Based on the concept that master equations can be derived from a microscopic 
description of the spin system interacting with external degrees of freedom, it is possible to gain 
insight into the underlying microscopic processes and to develop a systematic approach towards 
determining the specific form of reaction operator in concrete scenarios. 
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I. INTRODUCTION 

Spin chemistry is concerned with the interplay of chemical reaction kinetics and electron/nuclear spin 
dynamics. 1-3 The central, though not exclusive, paradigm of spin chemistry is the radical pair mechanism. 
It focuses on the chemical behavior of radical pairs, formed as reaction intermediates either as geminate 
radical pairs (RPs) in a correlated fashion from precursors of well-defined spin multiplicity, or in an uncor- 
rected fashion from the statistical encounter of free radicals (F-pairs). 4 The key assumption in explaining 
spin chemical effects is the postulate of spin conservation during elementary chemical processes such as 
electron transfer or covalent bond cleavage and formation, an assumption that leads to a dependence of the 
reactive behavior of a RP on the different spin substates. This notion implies that the relevant Hilbert space 
of the RP comprises a set of energetically close lying states differing essentially in electron/nuclear spin 
quantum numbers. (It should be noted that in case of the electron spin it is usually an effective spin, com- 
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prising the different spin-orbit mixed character of the various substates of the RP considered.) Insofar as 
the internal dynamics of this set of states can be described by an effective spin Hamiltonian, it is coherent. 
Hence, in a kinetic theory, the RP has to be represented by a density matrix p with the internal dynam- 
ics described by the commutator [H,p]. However, regarding the correct description of the spin-selective 
chemical kinetics, the question arises how to account for it within the Liouville formalism. Traditionally, 
in theoretical spin chemistry this has been done in an ad hoc fashion, but taking into account the necessary 
conditions for conserving the properties of a density matrix (positive, Hermitian). Very recently however, 
there has been a renewed debate about the correct way to treat the problem, with a number of different pro- 
posals ^ 10 ' 21 ' 22 , each starting with a different philosophy, using different formalisms and, more importantly, 
predicting different answers. Which of them is correct and which are wrong? Here we show a framework 
that allows to treat all presented approaches on the same footing. It demystifies the quantum measurement 
approach to spin chemistry, and provides a common basis, from which the various proposals can be ob- 
tained. Thereby, we can identify valid approaches, and we provide a systematic way of constructing new 
ones. The decision about which of the particular instances applies in a concrete scenario cannot be decided 
a-priori but depends on the precise physical situation. In more detail, the situation is as follows. 

For pedagogical reasons and to keep the analysis transparent, we will restrict ourselves to the case of a 
RP with exclusive singlet reactivity and, if concrete representations are given, we shall consider a simple 
two-level model (To, S) of the set of spin substates. The standard form of the reaction part of the stochastic 
Liouville equation (SLE) of the RP, as justified by the formal criterion that the density matrix must always 
remain positive definite, was given by Haberkorn 5 as 

\P\™ = -\UQsP+PQs) (i) 

where Q$ is the projection operator on the singlet RP state: 

Qs = \S)(S\ (2) 

and ks is the specific rate constant of the singlet RP reaction. 

Triggered by a paper of Kominis 6 who first proposed to treat the spin selective reaction of RPs in the 
framework of quantum measurement theory, there has been a renewed debate on the physically correct form 
of the reaction operator. Although the details of the formal derivation given by Kominis have been debated, 
the result 

[p]™ = -\h{QsP + pQs) + hQ s pQs (3) 

did arouse some interest. The difference between CO and © is that the latter conserves the trace of the 
density matrix. It describes a mere singlet/triplet dephasing due to the coupling of the RP states to the 
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reaction channel but does not account for the decay of the RP population, which Kominis took into account 
by a further ad hoc assumption of a rate of kgTr{QspQs} for the decay of RP population which, however, 
would not change the normalized density matrix of the RP 

Jones and Hore, 7 taking up the quantum measurement approach in a formally consistent way, arrived at 
the result: 

[p] rxn = -k s p + IcsQtPQt (4) 
= -\h{QsP + pQs) - \k s {QsP + pQs)+hQspQs- 

As indicated in the second line, using a resolution of the identity in the singlet/triplet subspace, Q$ + Qt = 1 
with Qj = \T)(T\ being the projector on the triplet RP state, this expression comprises the Haberkorn type 
reaction operator (first term) plus, in the second and third term the dephasing operator as it occurred in (|3), 
which is of Lindblad form. 

In follow-up work by Il'ichov and Anishchik, 8 Ivanov et al.j2 and Shushin 10 , it was attempted to 
derive appropriate forms of the reaction operator by starting out with explicit treatments of the interaction 
of the RP spin states with the environment, including the systems of intramolecular nuclear vibrations. All 
of these treatments essentially arrive at the Haberkorn reaction term and question the appropriateness of the 
quantum measurement-based view of the RP reaction. 

In this work, we start from the same premise that a master equation description for the state evolution 
originates from microscopic dynamics of the considered quantum system interacting with additional, usu- 
ally uncontrolled degrees of freedom named the environment. Although these degrees of freedom may allow 
for an approximate classical description, as a source of noise for example, a first-principle treatment should 
include them as quantum degrees of freedom. Therefore, the evolution of the composite system, i.e. system 
and environment, is coherent and described by a commutator that contains the system-environment interac- 
tion. When restricting the formal description to the system only, i.e. treating the open quantum system, this 
interaction results in an effective, in general non-unitary evolution that needs to be accommodated within 
the Liouville formalism, for example. From this perspective it is possible to identify classes of system- 
environment interactions that lead to the specific proposed reaction operators, on one hand. On the other 
hand, given that a master equation is an accurate description of the real reaction dynamics, this approach 
opens up new possibilities to identify the form of the reaction operator in concrete scenarios. 

Within the chosen framework, i.e. including the environmental degrees of freedom into the description, 
it is furthermore possible to phrase the open system dynamics in terms of measurements. Thereby we 
can embed the kinematics due to "quantum measurements" into a broader context, in which the quantum 
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system of interest interacts with environmental degrees of freedom. The effective dynamics of the system 
can always be interpreted as measurements by the environment. Along the lines that "any linear process 
affecting an open system A can be viewed as an unread generalized measurement, i.e. as resulting from an 
entanglement of A with an environment simulator B, followed by an unread projective measurement in B" 
as has recently been stated by Haroche and Raimond 11 and as it is known from measurement theory and 
open system dynamics—^—, we will show in this work that a systematic view on the various approaches to 
the reaction term in spin chemistry can be obtained by taking into account how the environment is affected 
through the interaction with the RP. 

We start out with defining dynamical maps that characterize in a very general way how the interaction 
may affect system and environment. From these maps we derive pertinent Lindblad operators in the Hilbert 
space comprising the spin states of the RP and the reaction product. Reducing these equations to the 
subspace of the RP leads to the various forms of master equations discussed in spin chemistry. Finally, we 
indicate possible experimental avenues for the characterization of the reaction operator. 

II. DYNAMICAL MAPS OF SYSTEM-ENVIRONMENT INTERACTION 

Before considering the dynamics, it is necessary to fix the concepts and to specify which degrees of free- 
dom are considered. The quantum system of interest are the spin degrees of freedom of the radical pair (RP). 
The spins are usually the effective spins of two unpaired electrons, each being localized on its respective 
molecule of the RP. In order to capture the full dynamics, it is necessary to take into account the different 
situations of whether or not the charges have recombined or a covalent bond has been formed. We therefore 
describe the system Hilbert space by the spin state of the unrecombined radical pair, which we limit for 
pedagogical reasons to the singlet state and the triplet state with zero projection of total spin, denoted by \S) 
and \T), respectively. In addition, we include into the Hilbert space the recombined state/reaction product 
state denoted by \P), which is also a singlet state but with a different electron configuration. (Formally, we 
thus consider a system of two degrees of freedom: electron spin and electron configuration. But in practice 
the reaction product will comprise an entire manifold of states orthogonal to the radical pair states \S) and 
\T). Since we do not intend to distinguish between different singlet reaction products, we can effectively 
lump together the entire reaction product manifold under the label \P).) 

The environment degrees of freedom shall remain abstract for the moment, since it is not necessary to 
restrict the treatment to specific examples. For concreteness, the reader may, however, think of the state of 
the electromagnetic field, which may be in the vacuum state or with a photon emitted due to fluorescence in 
the course of recombination, or the relative position and orientation of the radical pair molecules carrying 
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the unpaired electron spin, and the configuration of their solvent shells. 

Let us first consider the interaction between the system and the environment for some time interval 8t 
under the interaction Hamiltonian H mt . The quantum state evolution of system and environment can be 
captured completely by a map from all possible initial states to their respective final states. It is, however, 
sufficient to only consider all products of basis states of system and environment, thereby simply represent- 
ing the matrix elements of the time-evolution operator U (St) = exp(— iH mt 8t). The basis states of the two 
RP spins are the singlet state \S) S , and the triplet state with zero projection of spin, denoted by \T) S . In ad- 
dition to the basis states of the system, \S) S , \T) S , and \P) S , the environment, to which the system is coupled, 
is assumed to start from some initial state |0),,. For clarity, system and environment states are labeled by 
subscript s and e, respectively. 



A. Spontaneous decay 

The first example of a system-environment interaction is a straightforward realization of a recombination 
dynamics with a minimum number of terms, and we will recognize this to be a simple, spontaneous decay 
process. It is described by the dynamical mapping: 

\S) s \0)e^a\P)s\°)e + b\S) s \0) e , (5) 
\T) s \0) e ^ \T),\0) e . (6) 

The first part in © captures the fact that with probability \a\ an initial singlet spin state evolves into some 
reaction product state \P) S with the environment evolving into some corresponding state \a) e , e.g. with a 
photon emitted due to fluorescence. Note, however, that the environment state \o) e might be some effective 
state of many degrees of freedom comprising the environment. With probability \b\ = 1 — \a\ both remain 
unaffected. The triplet does not interact with the environment and hence both do not evolve under the 
interaction Hamiltonian. The probability amplitudes a and b depend on the coupling strength and the time 
duration 8t of the system-environment interaction because they are matrix elements of the time-evolution 
operator: 

a = e {a\ s {P\U(8t) \S) s \0) e , (7) 
b = e (0\ s (S\U(8t)\S) s \0) e . (8) 

This form of open system dynamics is generated by an interaction Hamiltonian of the generic form 

Hint =g(\P)s(S\ \a) e {0\ + \S) S (P\ \0) e {a\) , (9) 



6 



where g denotes the coupling strength. In the present work, we will only focus on the generic structure 
of the interaction Hamiltonian. For the discussion of specific examples that are relevant in the radical pair 
model see&~— . 

As a preparation for what follows, we need an estimation of b as given in ©. By expanding U (5?) to 
second order, and using the Hamiltonian d9j, we find that 

b « 1 -\ e (0\,(S\ Hi, \S) s \0) e 8t 2 = 1-e, (10) 

the second term of which we denote by e. Within this approximation e is a small positive quantity, and 
hence b is real. 

According to the dynamical map d5j and d6j, an arbitrary pure state that starts in the singlet/triplet 
subspace, i.e. with \cs\ 2 + \ct\ 2 = 1, would evolve as follows: 

(c s |S) 4 + cr|r),)|0) e -> (bc s \S) s + c T \T)^\0) e + ac s \P) s \o) e . (11) 

It is the properties of the combined state of system and environment that is important for the understanding 
of the intrinsic dynamics of the system. According to the axioms of quantum theory, a measurement that 
provides information of whether or not the reaction (with product state \P) S ) occurred in the system, yields 
a positive outcome with probability \acs\ and result in the measured state \P) s \o) e of system and environ- 
ment. A negative measurement outcome would simply return the (appropriately renormalized) state in the 
singlet/triplet subspace with the environment in |0) e . Since here a system state in the singlet/triplet manifold 
only appears together with environment state |0) e , and the product state only with \o) e , a measurement of 
the environment state instead of the system state would yield the same results. 

When only the system degrees of freedom are considered, it is necessary to trace over the environment 
degrees of freedom. Tracing out environmental degrees of freedom is formally equivalent to the already 
mentioned "unread projective measurement" interpretation. That is, a projective (von Neumann type) mea- 
surement of the environment state, e.g. of state |0) e or \a) e , also collapses the system into its respective 
state with corresponding probabilities, as given according to (fTTT) . for example. Ignoring the measurement 
outcomes is described by a mixed state composed of the individual states of the measurement outcomes, 
weighted by their respective probabilities of occurrence. The same result is obtained by the partial trace 
over the environment. A priori, this description applies to the quantum state of a single radical pair, i.e. the 
proper state of a single radical pair (molecule) is given by a mixed density matrix. In the case of ensembles 
of many radical pairs, measurement probabilities for the single radical pair translate into frequencies of 
occurrence when the same measurement is carried out on all members of the ensemble, and provided that 
there are no correlations between different radical pairs of the ensemble, and that the state of all radical 
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pairs are identical copies of each other. Only if the latter assumptions hold, the density matrix of the single 
radical pair provides a complete description of the ensemble. 

In the present context, it is of primary interest how the state evolves in the singlet/triplet subspace, and 
thus in addition to tracing out the environment, one projects the system onto this singlet/triplet subspace. 
The projection will in general yield a state p in the singlet/triplet subspace that is not normalized. The 
procedure of normalizing this state yields the (normalized) reduced state of the system that is produced by 
a measurement of whether or not a reaction product \P) e occurred, in the case of a negative outcome. The 
fraction of the state that started in an arbitrary initial pure state and has not yet undergone a recombination 
reaction is therefore described by the following mapping: 

|2 * \ , / 1 7 i2 i i2 



p(o)=l N csc \ hrr-nH N b "1 1 ='<*>■ 

c*c T \c T \ ) l-| flt? sl \b*c* s c T \c T \ 



Note that, as an obvious consequence of (fTTT) . a factor b appears only in combination with factors cs ■ 

An arbitrary, possibly mixed initial state of the electronic spin system will therefore obey the nonlinear, 
conditional evolution 

/ Pss Pst\ 1 / \b\ 2 pss bp S T \ 

P(°)= ^7-7^ \=P{8t), (13) 

\Pts PttJ 1_ M Pss \b*p TS Ptt J 

where the final state is normalized by the probability (1 — \a\ 2 pss) of not having reacted. 

In case of a single RP, this evolution would correspond to the evolution of the RP under the condition that 
it has not yet recombined. In case of uncorrelated ensembles of RPs, this translates into the subensemble 
of molecules that are still present (have not yet recombined) and whose absolute number scales with the 
normalization factor. 



B. Quantum measurement-induced pure dephasing 

A scenario different from spontaneous decay dynamics is given by a system-environment interaction 
without a recombination of the RP spins but where information regarding the singlet/triplet character of the 
RP spin is left in the environment, i.e. the spin character is effectively measured by the environment. In 
such a measurement interaction, the environment changes its state, depending on the RP spin character in a 
different way, which is captured by the following dynamical mapping for a time step 8t: 

\S) s \0) e ^a s \S) s \o) e + bs\S) s \0) e , (14) 
\T) s \0) e ^a T \T) s \r) e + b T \T) s \0} e . (15) 
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This dynamics is generated by an interaction Hamiltonian of the form 



flint = 8s\S)s(S\ ® (|«T) e (0| + |0) e (ff|) +g T \T) s (T\ ® (|T) e (0| + |0) e <T|) 



(16) 



where g$ and gj are the respective interaction strengths for the environment coupling to the singlet/triplet 
character. 

In contrast to the decay process in © and ©, the present mapping contains the additional term 
aT\T) s \x) e in (fT5T > for the triplet, and there is no state \P) S . Although the singlet and triplet components 
of the system do not seem to be affected, as the spin part of the composite final state is again a singlet 
or triplet with states \S) s (as\o) e + bs\0) e ) and \T) s [aT\t) e + £r|0) e ), respectively, the environment does 
change conditioned on the spin character and thereby also affects the spin state. With probabilities \as\ and 
\aT\ the environment senses the presence of the singlet and triplet, respectively, i.e. in particular the triplet 
character, which differs from the situation of the spontaneous decay dynamics. 

Thus, for the present interaction, the environment effectively performs a measurement of the spin be- 
cause it evolves conditioned on the spin character and thereby effectively senses the presence of a singlet 
or triplet. A measurement of the environment state would yield the respective state of the spin system, i.e. 
\S) S or \T) S for states \a) e and \x) e of the environment, respectively, or an undetermined state for |0) e . Ac- 
cordingly, the environment states \a) e and \x) e constitute the measurement outcomes of this singlet/triplet 
measurement, and thus play the role of pointers of a measurement device that indicates the presence of a 
singlet or triplet. Given a singlet or triplet state, this measurement produces a conclusive outcome with 
probabilities \as\ and |ar| , respectively. With the remaining probability \bs\ and \ot\ , the singlet or 
triplet state, respectively, is not detected, and the environment (measurement device) remains in its initial 
(neutral) state |0) e , which thus represents an inconclusive measurement outcome. 

Since during this dynamics, the singlet fraction in the subspace of the unrecombined radical pair does not 
change, - it is not taken out in the form of reaction products, and moved to the \P) S subspace - the singlet 
population is not affected by this dynamics. To simplify the following discussion, let us assume equal 
measurement probabilities ("detector efficiencies") for singlet and triplet. That is, we assume bs = bj = b 
and thus a$ = aj = a, and can therefore drop the index. When restricting our consideration to the spin 
system, as done for the decay scenario by means of tracing out the environment, the resulting spin state is a 
probabilistic mixture of the detected singlet/triplet fractions |S)j(S| and \T) S (T\, each with probability \a\ , 
and the undetected (coherent superposition) spin state that occurs with probability \b\ = 1 — \a\ 2 . Therefore, 
only the singlet/triplet coherences are destroyed during the evolution: 




(17) 
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When compared to the decay scenario (fT3T >. the state automatically remains normalized because it com- 

2 

pletely remains in the singlet/triplet manifold and, as a formal consequence, the factor \b\ in the singlet 
population disappears because here it is accompanied by the contribution \a\ of the detected singlet that is 
not taken to the reaction product state. 

This quantum measurement dynamics results in pure dephasing, but it does not capture the transfor- 
mation of the singlet into a reaction product. As such, pure dephasing might well appear as an additional 
process that takes place in the radical pair model and happens in parallel, but it cannot represent the recom- 
bination kinematics^. 

C. Quantum measurement-induced recombination 

In contrast to the spontaneous decay dynamics and pure dephasing treated in the previous sections, 
an interaction between the RP spin system and the environment is conceivable that involves two steps. 
First, the RP spin character is measured by the environment, and second, dependent on the measurement 
outcome, the detected singlet fraction is transformed into a reaction product. This can be accounted for 
by combining the dynamical maps of pure dephasing and spontaneous decay into a quantum measurement- 
induced recombination map: 

\S) s \0) e ^a s \P) s \a) e + b s \S} s \0) e , (18) 
\T) s \0) e ^a T \T) s \z) e + b T \T) s \0) e . (19) 

Here, we assumed that upon detection, the singlet is transformed with unit efficiency into reaction products 
and hence there is no term of the form |S). S . |<r) e . We emphasize that the fact, that the singlet fraction 
is transported into a reaction product state \P) S , after the environment has sensed the singlet, is rather 
irrelevant to the result that information about the singlet character is left in the environment in the form 
of the contribution \o) e . Note, that although we use the same notation for the states of the environment 
to indicate their correlation to the radical pair's spin state, here the character of the environmental state 
differs drastically from the case in ([Pfl i since product formation implies a change in the electronic orbital 
degree of freedom. Often during the product formation energy is released. One can thus generally expect 
that the changes in the environment due to the present type of measurement involve larger energy scales or 
reorganization processes than the dephasing treated in the previous section. 

To keep the following discussion transparent, let us again simplify the present scenario by assuming 
equal measurement probabilities (efficiencies) for singlet and triplet. That is, we assume b$ = bj = b and 
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thus a$ = aj = a. Then, an arbitrary pure initial spin state evolves according to 

'c s \s) s +CT\T) s )\o) e ^a(cs\P) s \a) e +c T \T) s \z) e )+b(cs\s) s +c T \T) s )\o) e . (20) 



2 

Thus, with probability \a\ , singlet and triplet are detected by the environment, which changes its state 
correspondingly, but in addition, the detected singlet fraction is transformed into the reaction product state 
\P) S . Again, with remaining probability \b\ 2 , initial spin state and environment remain unaffected. 

Please note, that in the (generic) case that the measurement probabilities are not equal, the spin state 
changes even if the state is not detected by the environment (because bs ^ bj cannot be factored out as is 
done with b in the equation above). Therefore, the state evolution under unsuccessful measurements by the 
environment will also have physically observable consequences ! 

By tracing out the environmental degrees of freedom and projecting to the singlet/triplet manifold, the 
evolution of the density matrix of the initial spin state in the singlet/triplet subspace is described as 



w 




K i 1 A 


\C*CT 


\ c t\ 2 ) 


' l-\ac s \ \ 



,2 I ,2 i,i2 * 

| \c s \ \b\ c s c*j 



(21) 
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where the final state appears only with probability 1 — \acs\ of the state not having decayed. This is the 
normalization factor. Note that in contrast to the spontaneous decay scenario in (fTTb . where a factor b was 
associated with every occurrence of cj, now a factor b is also associated with every occurrence of cj. In 
the triplet population, the term \a\ + |&| 2 = 1 appears as a pref actor. An arbitrary mixed initial state of the 
electronic spin system thus evolves according to 

, v , , , 1 / \b\ 2 pSS HVsA . 

P(0)=( , \=P(St), (22) 

1 - \a\ Pss \\b\ prs Ptt J 

2 

where, again, the final state is normalized appropriately, but here by a factor 1 — \a\ pss- This equation cap- 
tures the conditional evolution of the uncorrected subensemble of radical pairs that have not yet recombined 
during the time step 8t. 

A comparison between the spontaneous decay scenario ([T3T l and the quantum measurement-induced 
recombination scenario (1221 shows that, at the level of the spin system, the only difference is that the 
coherences in the latter scenario carry the factor \b\ instead of b or its complex conjugate. However, from 
the viewpoint of what happens to the environment, the underlying physics differs drastically. The mapping 
for the present reaction kinetics in (TT8T i and (fT9l l is generated by an interaction Hamiltonian of the form 

Hint = gs(\P)s{$\ ® \o) e (0\ + \S) S (P\ |0> e (ff|) +g T \T),{T\ 8) (|T) e (0| + |0) e (T|) . (23) 
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Compared to the interaction Hamiltonian of the spontaneous decay, the difference in comparison to the 
present case is what happens to/due to the triplet character of the state! On one hand, the interaction of 
the triplet with the environment induces additional dephasing in the reduced state of the spin system, on 
the other hand, the environment state is changed due to the presence of the triplet. Since here the reaction 
products originate only from the singlet character of the radical pair, one could think that the reduced 
dynamics of the spin systems is only affected by the singlet character of the state. However, in addition, 
the open-system dynamics of the triplet, which does not recombine, influences the state dynamics. From 
the viewpoint of the mappings and the interaction Hamiltonian, the present quantum measurement-induced 
recombination dynamics can thus be seen as the combination of a spontaneous decay and a pure dephasing. 

Every conceivable interaction between the system and an environment can be represented in terms 
of certain general types of evolution mappings (so-called completely positive maps, where one assumes 
that the initial state is given as a product state between the system and the environment)^, out of which 
we have chosen three relevant examples. In addition to the mappings described above, alternative and 
more complicated recombination dynamics can be constructed, which may capture more effects and would 
include more states of the environment. Ultimately, it remains to the experiments to distinguish, which 
physical degrees of freedom are relevant and which roles they play for the recombination dynamics. 

An important lesson that can be taken from the presented examples is that an alternative way of distin- 
guishing spontaneous decay reaction kinetics from quantum measurement-induced reaction kinetics, is not 
only by the rate at which coherences of the RP spin state evolve, but by the traces that have been left in the 
environment. 



From the interaction between system and environment that is described in the previous section by a 
unitary mapping for a single time step 8t, one can derive, using a coarse graining in time and additional 
assumptions, master equations of Lindblad form. 1516 Such master equations can also be obtained within the 
original Bloch-Redfield approach to describe magnetic relaxations^^. The general form of a Lindblad- 
type master equation is given by 



for the incoherent part of the evolution. The operators L, are called Lindblad or jump operators, and are 
non-hermitian in general. The k, are rates and are positive. (Since we focus on the system-environment 
interaction only, the coherent part of evolution due to the system Hamiltonian — i[H s ,p] is omitted, which 



III. MASTER EQUATIONS 




(24) 
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effectively amounts to using an interaction picture or simply setting H s = 0.) This form of evolution equation 
guarantees that the density matrix remains normalized and positive. The latter means that the probabilities 
(eigenvalues of the density matrix) are non-negative. The master equation is time-local, i.e. p(t) depends 
only on p(t), and therefore cannot incorporate information of the past into the evolution. In a rigorous 
derivation, this memoryless evolution results from the Markov approximation, where a time integration 
over environment correlation functions for earlier times is replaced by the properties of the environment at 
the current time. 

If one wants to attribute a functional interpretation to the individual terms in the Lindblad master equa- 
tion, the first term in (l24l i can be interpreted as the state change due to a "successful" interaction with the 
environment, whereas the second term corresponds to the state change due to an "unsuccessful" interaction 
event, i.e. the part of the state that remains unchanged due to the system-environment interaction but is 
renormalized by the corresponding probability and thereby changed. In the present context, the first term 
will be associated to the recombination kinematics that transforms the radical pair to the reaction product, 
and thereby takes population out of the singlet/triplet manifold. The second term leaves the populations 
of the radical pair spin characters invariant, but modifies the conditioned evolution, which thus exhibits a 
non-unitary component. 

The mappings introduced in the previous section capture the action of the time evolution operator U for 
a short time St. The expressions on the right hand side of these mappings thus represent the result of U (St) 
applied to the initial state, e.g. for a pure state: 

\V(5t)) se = U(5t)\V(0)) se (25) 
= exp(-iH int 8t)\xi/ ) s \0) e . (26) 

Each repeated application of the mapping transports the state forward in time by an additional time incre- 
ment 8t. If we assume that in every application of the mapping, the system interacts with a fresh initial 
state of the environment [0) e , we neglect memoiy effects of the environment due to previous interactions, 
i.e. we perform the Markov approximation. This assumption is often reasonable if there are sufficiently 
many environmental degrees of freedom available, and/or the internal dynamics of the environment is fast 
enough such that correlations decay quickly and the environment effectively returns to its previous state, 
e.g. a photon or heat released into the environment quickly propagates away and does not come back. 

In order to compare the dynamics on the level of master equations, rather than on the level of dynamical 
mappings, we will now discuss the form of the Lindblad master equations arising for the three paradigmatic 
cases of reaction kinematics presented in the previous section. 
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A. Spontaneous decay 

From the dynamical map of the spontaneous decay kinematics, we can derive the differential evolu- 
tion equation for the density matrix of the system. Instead of deriving this master equation directly from 
the system-environment interaction Hamiltonian, we employ the dynamical maps and provide a sketch of 
how to relate them to the master equations of Lindblad form. For detailed derivations of Lindblad master 
equations from system-environment interactions see Refs. UilfL, for example. 

We start with the unitary map applied to a general system state as given in (fTTT > and trace over environ- 
ment degrees of freedom: 

P (0) P (St) = \a\ 2 \p){s\p(o)\s)(p\ + (b\s)(s\ + |r)<r|) P (o)(V|s><s| + |r)(r|). (27) 

Since p(0) is assumed to have only components in the singlet/triplet manifold, expressions of the form 
\T)(T\p(0) can be replaced with p(0) — |S)(S|p(0) and we obtain, using the notation Q$ = \S)(S\ and Qj = 
\T)(T\ for the projectors: 

p(8t) = p(0) + \a\ 2 \P)(S\p(0)\S)(P\ + \b-l\ 2 Q s p(0)Qs + (b-\)Q s p(0) + (b* - \)p(0)Q s . (28) 

Since, for short times 8t, we can write b = 1 — £ with a small real quantity e (cf. equation (fTOl)). it follows 
from \a\ 2 + \b\ 2 = 1 that \a\ 2 « 2s in leading order. With these relations, and keeping only terms linear in e, 



equation (1281 ) is reduced to: 

Sp=p (St) - p (0) « \a\ 2 \P) (S\p(0) \S) (P\ + (b - 1) ( Q s p (0) + p (0)fis ) . (29) 



Comparing this with the general expression for the Lindblad master equation (I24I) . we see that there is only 
one Lindblad operator, with the following assignments: 

L=\P)(S\, (30) 
L^L = \S) (P\P) (S\ = \S)(S\ = Q s , and (31) 
\a\ 2 = kSt (32) 

which is consistent with \a\ 2 + \b\ =1 and hence b w 1 — \ \a\ for small 8t. The resulting Lindblad 
operator is also suggested by the interaction Hamiltonian © because it induces a "quantum jump", i.e. 
a change of state from the singlet spin component to the reaction product state with rate k, and thereby 
yields the standard master equation for spontaneous decay dynamics 1 -. For example, for the case of an 
electromagnetic field mode in a leaky cavity, L would simply be the corresponding photon annihilation 
operator, while for an excited atom that interacts with the electromagnetic vacuum, L would be the atomic 
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lowering operator. In the context of the radical pair mechanism, Gauger et al. have recently employed this 
form of a spontaneous decay to model the recombination dynamics. With the expression for L at hand, the 
master equation (l24l takes the form 



After performing a projection to the subspace of the singlet/triplet component, the first term vanishes, and 
the master equation for the reaction operator reads 



which is just the Haberkorn form of the recombination dynamics in (fl]). Since the first term of the Lind- 
blad equation is missing, it can no longer preserve the trace of the density matrix within the singlet/triplet 
manifold. In matrix form, the equation reads 



It shows that the coherences in the singlet/triplet basis decay at half the rate of the singlet population. A 
decay out of the triplet manifold can be included in an analogous way. 

At this point a short remark about coherences and their decay is in order. In the context of radical pairs, 
coherences are always considered in the basis of singlet and triplet, i.e. matrix elements of the kind (S\p\T). 
Note, that in other fields coherences in a quantum system are usually considered in a product basis of its 
constituents. The singlet/triplet basis states are eigenstates of a collective property of the coupled spin pair, 
whereas the product basis {|tt)> |t4-)> 1 4-T) ? 144-)} is represented by the outer products of basis states of the 
individual (uncoupled) spins. The loss of coherences (dephasing) in the singlet/triplet basis, for instance due 
to a singlet/triplet measurement, is thus not equivalent to a dephasing in the product basis, i.e. the decay of a 
matrix element of the kind (ti\p |4t)- The latter emerges due to local measurements of the individual spins, 
for example. Whereas coherence of the former type represents the presence of (electron) spin polarization in 
spin chemistry, i.e. a preponderance of | j) or \\) electron spin on the individual radical species, coherence 
of the latter type goes along with the concept of entanglement in quantum physics. 




(33) 




(34) 




(35) 



B. Quantum measurement-induced pure dephasing 



In a similar fashion as it was done in the case of spontaneous decay dynamics, one can construct the rela- 
tion between the unitary mapping and the master equation for the pure dephasing dynamics. For dephasing 
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in the singlet/triplet basis, a corresponding argument leads to the two Lindblad operators 

L S = \S)(S\ = Q S and L T = \T)(T\ = Q T , (36) 

i.e. the measurement operators of the singlet and triplet state. With respective rates k~s = \as\ 2 /St and 
kj = \aj\ I St, the Lindblad equation then reads 

[Plrxn = L ** [QiPQi -UQiP+PQi)) , (37) 

i=S,T V z / 

= fc+h)(&pQs-^(QsP+PQs)\ (38) 

where for the second step we used that Q$ + Qt = 1 in the singlet/triplet subspace. According to this master 
equation, the density matrix evolves as 




-\{ks + ~kr)psT 



(39) 



k Pts Ptt J \-2( k s + h)PTS 
Here, only the coherences decay with rate (ks +kr)/2. When compared to the dephasing master equa- 
tion ©, this rate amounts to 

~k s + ~k T =k s , (40) 

i.e. the effective rate appearing in ©. 

Regarding the additional ad hoc expressions for the formation of reaction products as introduced by 
Kominis in various ways^i, it is unclear whether a microscopic physical model with a system-environment 
interaction exists from which these effective, phenomenological master equations can be derived. 

C. Quantum measurement-induced recombination dynamics 

Finally, we construct the Lindblad-type master equation for a quantum measurement-induced recom- 
bination dynamics. From the infinitesimal change of state in the system, which we obtain from the map- 
ping (l22l ) in the full space, we arrive at the following equation when proceeding in a similar way as explained 
above: 



S P = P (St)-p(0)^\a s nP){S\ P (0)\S)(P\ + (b s -\){\S)(S\p(0)+ P (0)\S){S\) (41) 
+ \a T \ 2 \T)(T\p(0)\T){T\ + (br - l)(\T)(T\p(0)+p(0)\T){T\) . 
We can now identify two Lindblad operators, which are also suggested by the interaction Hamiltonian d23l : 

L S =\P)(S\ with \a s \ 2 = k s St and (42) 

L T = \T)(T\ = Q T with \a T \ 2 = k T St. (43) 
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After projecting again to the singlet/triplet subspace, the master equation no longer contains the term L$pL s 
and, by using L^Lg = Qs as in (I3TI ). it takes the form 

[AU = ~\h(QsP + pQs) +k T (qtPQt-\(QtP + PQt)^ (44) 
= - \h {QsP +PQs)+ hQspQs - \h {QsP +PQs) (45) 
= - ] ^ L (QsP+pQs)+k T QspQs. (46) 

For the special case when the singlet and triplet component are measured at the same rate (i.e. through 
a mechanism with the same efficiency), k$ = kj, the second line reduces to the Jones-Hore master equa- 
tion (01). The individual components of the density matrix now evolve according to 

d I Pss Pst\ _ I -ksPss -j(ks + k T )psr 
dt \Pts Ptt) \-^(ks + k T )pTS 

where in the case k$ = kj the coherences decay at the same rate as the singlet population, i.e. at twice the 
rate of the spontaneous decay in eq. (I35T ). 

Looking at the underlying physical processes, i.e. at the full space, we see again that the relevant differ- 
ence between the spontaneous decay scenario and the measurement-induced recombination is the dynamics 
that affects the triplet state. It relates to the information about the triplet character that is obtained by the 
environment through an interaction or "unread measurement". For a suitable choice of rates k$ and kj, one 
obtains an entire manifold of processes that may show various decay rates of the coherences. In the limit 
kj = 0, i.e. vanishing interaction of the triplet with an environment and thus no measurement of the triplet 
character, we obtain the Haberkorn recombination dynamics of a spontaneous decay. In the other limit of 
ks = 0, a pure dephasing dynamics emerges. Intermediate regimes with both k$ > and kj > yield re- 
combination kinematics of the kind similar to the quantum measurement-induced recombination, with the 
Jones-Hore kinematics as a special case for k$ = kj. 




IV. DISCUSSION 



Let us finally compare the three considered models for the reaction kinetics. It is obvious that, on the 
level of the system in the RP singlet/triplet subspace, three entirely different system-environment interac- 
tions create very similar time evolutions that only differ in the rate at which the coherences decay. Observing 
only the RP spin system in the singlet/triplet subspace may therefore not be sufficient to decide on the gen- 
erating process and the underlying physics of the recombination kinematics. For example, by combining 
several independent processes affecting the radical pair at the same time, their Lindblad master equations 
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are simply added. Therefore, in a process where, on one hand, a spontaneous decay from the singlet at rate 
k$ and, on the other hand, a pure dephasing ("measurement") process with rates k$ = ks/2 and kj = ks/2 
affect the spin system simultaneously, the evolution can be interpreted as a quantum measurement-induced 
recombination with rates k$ = kj as proposed by Jones and Hore. Indeed, the master equations for the 
radical pair would be identical. It is therefore not surprising that a similar phenomenology on the level of 
the RP spin systems can be obtained by even more complicated models as recently considered by Jones et 
alj22, which include a larger number of reaction stages than those considered here. 

In order to distinguish between dynamics due to a quantum measurement-induced recombination, and 
a suitable combination of spontaneous decay and dephasing, it is necessary to also consider the involved 
degrees of freedom of the environment. Including the environment degrees of freedom into the description 
allows us to identify the origin and influence of "measurements" during the recombination dynamics. The 
dynamics caused by the interaction of system and environment represents a situation in which the environ- 
ment effectively measures the (spin) state of the system. Since these "measurement results" are never read 
out, the single radical pair system per se is correctly described by a mixed state that evolves according to the 
master equation. We thereby complement previous treatments 8 - ~— focusing on the problem of identifying 
possible relevant interaction Hamiltonians, by embedding them into a framework that allows one to identify 
quantum measurements of the radical pair spin system. 

By including the interaction with a (quantum) environment from the beginning, the origin of the terms 
in the master equation becomes transparent, and avenues emerge that allow for distinguishing the various 
proposed reaction mechanisms by the traces that the state of the system leaves in the environment. Once 
the relevant environmental degrees of freedom are identified, their experimental examination would provide 
insight as to what information about the RP spin character has been left in the environment. The absence 
of an interaction, however, can only be stated with some uncertainty determined by the experimental pre- 
cision. Furthermore, an experimental measurement of the environment would disclose information on the 
system, and thereby project it into the respective system state. It would thus become possible to unravel the 
decoherence dynamics of the RP spin state, and observe the system dynamics as it stochastically evolves 
under the system-environment interaction. Such an observation realizes the dynamics in the form of indi- 
vidual stochastic (almost) pure state trajectories instead of the averaged mixed state description provided 
by the master equation. Similar processes are well-known in the field of quantum optics, where suitable ex- 
periments have been performed with systems of trapped ions 23 , micro-masers and cavity QED-setups^i 2 ^. 
A translation of these or similar experiments to single molecule (i.e. single radical pair) systems of spin 
chemistry (or quantum dot systems that might mimic radical pairs) represents an experimental challenge, 
but would yield an unprecedented level of access to spin chemical reactions on the quantum level. Possible 
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experimental handles to environment degrees of freedom may include single photon detection for fluores- 
cent radical pair recombination, the detection of dissipated heat by spectroscopy or thermal lens/thermal 
grid techniques, and the resolution of the intermolecular Brownian molecule motion. 

In the present treatment, we have restricted ourselves to a simple "two-level system" for the RP spin 
manifold consisting of the \S) and \Tq) spin states. This scenario is approximately realized for singlet- 
born radical pairs in the high-field limit, where the spin states \T±) are energetically separated due to a 
large Zeemann interaction. However, possible open system dynamics of a two-level system are limited to 
essentially T\- and 72-type processes. For the interesting low-field regime, however, the fully accessible 
four-dimensional spin space offers a richer dynamical phenomenology. 

V. CONCLUSION 

The presented framework gives a common basis, on which the various reaction operators proposed in 
the literature can be equally constructed and compared. By adopting a view of system and environment, 
it becomes natural to identify quantum measurements in the reaction processes, thereby demystifying the 
quantum measurement approach to spin chemistry. Within the presented framework, we exemplarily derived 
the reaction operators proposed in the literature, while providing a systematic way of constructing new ones. 
Although both the traditional as well as the "quantum measuremenf'-inspired reaction operators emerge as 
formally correct proposals, the ultimate decision about which of the reaction operators are the appropriate 
model in a given scenario will certainly depend on the precise physical situation. 

The open quantum system approach to the spin-dependent recombination kinetics of radical pairs pre- 
sented here provides a general framework for justifying and establishing appropriate reaction operators 
for the stochastic Liouville equations customary in spin chemistry. Proceeding from general evolution 
mappings, that specify the particular interactions between spin system and environment, through pertinent 
Lindblad equations in the complete system/environment space to the final stochastic Liouville equation of 
the reduced spin system, defines a transparent approach to the theoretical treatment of spin-dependent rad- 
ical pair kinetics. An important lesson that we can take from this discussion is the following: We see that 
the rate at which the coherences between singlet and triplet decay is no longer the primary criterion to de- 
cide whether a quantum measurement of the spin character occurred. The discussion of the correct master 
equation should therefore not focus on formal expressions of rates but rather include the traces of the spin 
character that have been left in the environment. It is the involvement of the environment that justifies the 
description and treatment of radical pair kinetics in terms of quantum measurement theory. The interpreta- 
tion of the recombination kinetics in terms of quantum measurements can thus be embedded in the general 
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context of open quantum systems. 
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